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Partially-projected Gutzwiller variational wavelunctions are used to describe the ground state ol 
disordered interacting systems of fermions. We compare several different variational ground states 
with the exact ground state for disordered one-dimensional chains, with the goal of determining a 
minimal set of variational parameters required to accurately describe the spatially-inhomogeneous 
charge densities and spin correlations. We find that, for weak and intermediate disorder, it is 
sufficient to include spatial variations of the charge densities in the product state alone, provided 
that screening of the disorder potential is accounted for. For strong disorder, this prescription is 
insufficient and it is necessary to include spatially inhomogeneous variational parameters as well. 

PACS numbers: 71.10.Fd,71.23.An,71.27.+a, 
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I. INTRODUCTION 

The problem of how to treat theoretically systems of 
interacting electrons in disordered materials extends over 
many decades. The most difficult problems to solve are 
ones in which the physical properties of interest depend 
on spatial correlations of individual realizations of the 
disorder potential. In this case, the majority of disorder- 
averaged approximations fail and more sophisticated an- 
alytical or numerical techniques are required. Anderson 
localization and Coulomb gap physics are two examples 
of such problems^. The problems are particularly chal- 
lenging in narrow band materials where conventional ap- 
proximations for the Coulomb interaction are poor.^ 

There has been particular interest in low- 
dimensional disordered systems in recent years. 
For example, many transition metal oxides, espe- 
cially the quasi-two-dimensional high temperature 
superconductorsr^i^i^iiiSiiSiiiS are intrinsically disordered 
by chemical doping, and have electronic properties 
that appear to be strongly susceptible to disorder. 
More recently, experiments have demonstrated the 
formation of a two-dimensional electron gas at the 
interface of an otherwise insulating La2Cu04/Sr2Ti04 
heterostructurejiiii^ii^ii^ and there are indications 
that the electronic states at the interface are strongly 
influenced by inhomogeneitiesj^^ A more fundamental 
question that has arisen in the past decade is whether 
interactions may drive a delocalization transition in thin 
metal films. 

Theoretically, we can investigate the relation between 
disorder and electron-electron interactions by studying 
the Anderson-Hubbard (AH) Hamiltonian. The AH 
model is the Hubbard modelii for a lattice with site en- 
ergies chosen from a random distribution 
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where i,j = 1, . . . , TV denote the sites of the lattice, (i, j) 
implies that i and j are nearest neighbours, c^.o- {fi-i,a-) is 



the destruction (number) operator for an electron at site 
i with spin cr, and the hopping energy is denoted by —t. 
The on-site energy at site i is given by Vi, and U is the 
electron repulsion of two electrons sharing the same site. 

The only known exact solution for the AH Hamilto- 
nian is in infinite-dimensionSfi^ and we thus have to 
use approximate methods in finite dimensions. The 
various approximate methods used previously include 
self-consistent Hartree-Fock theory;i2i2£iSii2^i^ dynami- 
cal mean field theories r^'^^iS^i^i^^iSS'^i^ and exact nu- 
merical calculations for small clusters i^i^^i^i^i Of par- 
ticular relevance to the current work is a recent varia- 
tional wavefunction approach by Pezzoli et al^ which 
we discuss below. 

Here, we consider a variational wavefunction approach, 
based on a modification of the Gutzwiller wave func- 
tion (GWF)^i22i^ to include spatial inhomogeneity. 
The GWF is the simplest variational wave function for 
Hubbard-type Hamiltonians and is given by 



I^'gwf) = - (1 - ff)n-»T^u)lV'ps) 



(2) 



where IV'ps) is a reference product (i.e. Slater- 
determinant) state. This function was originally intro- 
duced to study the correlations of the ground state of 
the Hubbard Hamiltonian. The variational parameter 
< g < 1 incorporates the effect of the Hubbard repul- 
sion between electrons of opposite spins on the same site, 
and is obtained by minimizing the energy functional 
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We note that g = I corresponds to an unprojected wave- 
function, while g = corresponds to a fully projected 
wavefunction in which there are no doubly-occupied sites. 
Physically, the g = projection captures an essential 
feature of the large-[/ Hubbard model: it generates lo- 
cal paramagnetic moments without breaking the spin- 
rotational invariance of the lattice. 

The GWF is hard to treat analytically, even for 
the disorder-free Hubbard Hamiltonian where analyti- 
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cal results have only been obtained in one^i^i^i^ and 
infinite'*'^ dimensions. The variational GWF is therefore 
primarily a numerical method.— "^"^i^ Our goal is to 
investigate the quality of a number of simple variational 
wavefunctions by comparison of the variational ground 
states with exact diagonalization calculations on small 
clusters. 

The GWF has certain well-known limitationSf2§, it fails 
in the disorder-free case, for example, to describe the 
Mott transition at half-filling in finite dimensions. In gen- 
eral, the GWF can be improved by the addition of nonlo- 
cal Jastrow factors that describe the long-range interac- 
tions between charge-density excitations.— Recently, 
Pezzoli et. al^ used variational methods to study the 
AH Hamiltonian in two dimensions. Their particular fo- 
cus was on states near the Mott transition, which necessi- 
tated the inclusion of the Jastrow factors. Furthermore, 
they included spatial inhomogeneity in both the prod- 
uct state iV'ps) and in their local variational factors (e.g. 
9 ~^ 9i)- These extra degrees of freedom resulted in a 
large number of variational parameters (several hundred 
for a typical lattice) and required sophisticated minimiza- 
tion schemes for the energy functional4^ It is therefore 
worth asking under what circumstances we can safely re- 
duce the number of variational parameters and still ob- 
tain reasonable results for the variational wavefunction. 

In the current work, we examine trial wavefunctions 
with a small number of variational parameters. Spa- 
tial inhomogeneity is incorporated in the product state 
wavefunctions IV'ps), however, the variational parameters 
are taken to be spatially homogeneous. Our main result 
is that for weak and intermediate disorder, this simple 
wavefunction gives a surprisingly good description of the 
ground state away from the Mott transition. It is only 
for strong disorder that our simple ansatz for I^'gwf) 
breaks down. 

The paper is organized as follows: In Sec. |TT1 we de- 
scribe the variational wavefunction in detail, and intro- 
duce a number of candidate states for \i'ps)- In Sec. IIIIl 
we asess the quality of the different trial wavefunctions 
by comparing with both exact diagonalization and un- 
restricted Hartree-Fock (UHF) calculations. The lat- 
ter comparison is motivated by the fact that UHF has 
been, until recently, the standard numerical technique 
for studying disordered systems. Furthermore, a recent 
work^" has suggested that UHF actually provides good 
quantitative results for some local physical quantities in 
disordered systems (namely that disorder improves the 
quality of the UHF approximation). We discuss the 
strengths and limitations of the simple GWF in Sec. llVi 
and conclusions are given in Sec. IVl 



want to keep the simple form of the projection operator 
in Eq. ([2]), which means that spatial inhomogeneity is 
introduced entirely through the product state. It is for 
this reason that we consider a variety of product states, 
with the goal of determining how large an effect the initial 
choice for IV'ps) makes on the final |5'gwf)- 



A. The Product States 

Assume we are seeking the A'^e-electron variational 
ground state for an A'^-site lattice with a particular re- 
alization of the disorder potential. We then take IV'ps) 
to be the A"e-electron ground state of an A^-site bilinear 
(effectively noninteracting) Hamiltonian TYbi whose dis- 
order potential is determined by the site energies of Ti,. 
We examine three possibilities: (i) Tihi is the Anderson 
Hamiltonian, obtained by setting [/ = in Eq. ([1]) and 
renormalizing the site energies with a variational screen- 
ing parameter, (ii) Tibi is the self-consistent Hartree-Fock 
(HF) decomposition of 7i, (iii) Tibi is the HF decomposi- 
tion of Ti with screened site energies. 

Since TLhi is bilinear, IV'ps) can be represented by 
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where a labels the lowest energy single- particle eigen- 
states of Hbi ordered from lowest [a — 1) to highest 
{a = A'e) energy, and 7^'s are creation operators for 
these eigenstates. The 7j,'s can be related to the real- 
space creation operators, c|, by: 
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where Gi^a are elements of the unitary matrix Q_ 

that diagonalizes Tihi, I labels the spin and site with 
1 < / < 2A. Then we have 

2N 

IV'ps) = E Gj,,^...Gj,^,nA,---c\JO), (6) 

/i,...,/n^ = 1 

which can be rewritten, using the anti-commutation re- 
lations of the creation operators, as 

l^ps)= E D[GU)]ci...cljO), (7) 

/l<...</iVe 

where D[G_{I_)] is a Slater determinant 



II. METHOD: GUTZWILLER VARIATIONAL 
APPROACH 

In this section, we introduce a number of variational 
states I^'gwf) that will be used for calculations. We 
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The vector 7 = ( Ji , . . . , 1^^ ) labels different configura- 
tions of the electrons. For a given /, gives the 

weight of that configuration of electrons in the product 
state. Using Eq. we can write 

|*gwf) = (9) 
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DFSGW 


Uniui) = 


VI = V.ls 


PMGW 


Uui^ = Uuii 


VI = V, 


PMGW(ff,e) 




V' = V/e 



TABLE I: Summary of the different bilinear Hamiltonians 
Tihi used to generate \ipps)- 



N 

E DML)] 11(1 - (1 - 5)«.T^U)4, ■ • ■ 4, |0)- 

/l<...</iVe j = l 

In general, I^'gwf) contains a large number of terms 
and must be calculated approximately, for example using 
variational Monte Carlo methods. In the current work, 
we restrict ourselves to small clusters where I^'gwf) can 
be evaluated exactly. 

Our first choice of Tibi is the non-interacting disordered 
Hamiltonian where the disorder potential is screened. 
The on-site energies are reduced. 



hare disorder potential ^ 

screened disorder potential 

where the screening factor e is a variational parameter. 
Both variational parameters e and g are determined so 
as to minimize the total GWF energy i?GWF(£,5)- We 
refer to this variational wavefunction as the disordered 
fermi sea GWF (DFSGW). 

Our remaining choices of Hbi are based on the HF 
decomposition of Ti. The four-fermion interaction term 
in Eq. ([T]) can be written as 

Uni]'n^l = U{niif + Sniif){nti + Sh,i) (10) 

where Uia- = {fiia), and dhia = hi^ — Uia-i^ The HF 
Hamiltonian, using Eqs. ^ and (fTO)) . can be written as 

N 

TLhf ^Hu-U y^^n^inii, (11) 
1=1 

where the second term on the right-hand side is just a 
constant, and 

Hbi = J2 + Un^-a)n^a - t ^ 4cj,, (12) 

i <y (ij)a 

where determined self-consistently. 

Our second choice of Tibi is the paramagnetic HF 
(PMHF) Hamihonian, i.e. Eq. with = n^. 

Here g is the only variational parameter, and we refer to 
the variational wavefunction for this choice as PMGW. 

Our third choice of Tibi is the PMHF Hamiltonian with 
a screened disorder potential. For this case, the on-site 
energies of Eq. are reduced to ¥( — Vi/e, and again, 
both variational parameters e and g are determined so as 
to minimize the total GWF energy. These calculations 



are more computationally demanding than the previous 
two because nio- must be determined self-consistently for 
each £. We refer to the variational wavefunction for this 
choice as PMGW((7,£). The different Tibi are summa- 
rized in Table |T1 

As a simple illustration, in Appendix|A]we compare the 
different product states for a disordered two-site system. 
As in the disorder-free case^^ the sample GWFs span 
the Hilbert space of a the two-site case and yield exact 
results. 



III. RESULTS: COMPARISONS OF EXACT 
AND VARIATIONAL QUANTITIES 

In this section we show results for PMGW, 
PMGW(5,£), and DFSGW variational states for small 
clusters where variational calculations may be compared 
to exact diagonalization calculations. Motivated by the 
earlier work^S demonstrating that the UHF approxima- 
tion (in which the self-consistently determined fiicr may 
depend on a) provides accurate charge densities in dis- 
ordered systems, we also consider UHF as a benchmark 
for the various GWF calculations. 

The quantities that we have calculated are as follows: 

• We have calculated the absolute difference of the 
exact and the variational energies per site defined 

by 

(5£;var = ^|S™-S<='^|. (13) 

• We have evaluated the magnitude of the overlap 
between the exact and variational wave functions, 
given by 

K^cxl^-var)!. (14) 

• We have calculated the local charge densities ac- 
cording to 

= (^l^n,,,!^'), (15) 

cr 

where \^) represents the exact or variational wave 
function. 

• We have calculated the local spin correlations for 
near-neighbour sites in one dimension, 

(S,-S,+i) = (^IS.-S.+il^'). (16) 
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FIG. 1: (Colour online) Variational parameters for the DF- 
SGW and PMGW approximations for the six-site cluster for 
weak {W/t = 4/3), intermediate {W/t = 4), and strong disor- 
der {W/t = 12). The DFSGW (a) screening factor e™^i„, and 
(b) projection gmin are shown as functions of U. The PMGW 
projection gmin is shown in (c). 



• We have calculated the average absolute difference 
of charge densities given by 

1 ^ 

1=1 

• We have calculated the average absolute difference 
of local spin correlations defined by 

(<5S-S')^-^ |(S..S,+i)™^-(S,.S.+i)'"^|. (18) 

i=l 

A. Six-Site Cluster 

We have examined several complexions of disorder for 
the 6x1 cluster with periodic boundary conditions at 
half-filling {Ne = 6), and focus on one representative con- 
figuration. The disorder potential for the configuration 
is 

V, = 0, -0.18W, +0.5W, +0.121^, -0.51^, +0.3W, 

where W is the strength of disorder. The bandwidth 
D = At for the noninteracting ordered cluster sets the 
scale for the disorder potential, and we have examined 
weak {W/t = 4/3), intermediate {W/t = 4) and strong 
{W/t = 12) disorder. Our results are shown in Figs. [T]to 

El 

We show first, in Fig. [U the variational parameters 
5mm and Smin that minimize the total energy for the 
DFSGW and PMGW approximations. The PMGW(5, e) 
approximation is numerically equivalent to PMGW (i.e. 



£min = 1) for weak and intermediate disorder, and is 
therefore not shown. 

For the DFSGW approximation, Smin is an increas- 
ing function of U , indicating that interactions effectively 
screen the disorder potential, making the charge distribu- 
tion more homogeneous than in the noninteracting case. 
At weak and intermediate disorder, Smin is only weakly 
affected by the disorder potential; however, for strong 
disorder there is a significant reduction in screening. For 
all W, Erain U when U > W. For weak disorder, 
<7min is quantitatively like that of the ordered case^^ 
while for strong disorder the wavefunction is unprojected 
(ffmin = 1) for U < 3.4t. In all cases, the projection is 
substantial when U ^ W,D. 

For the PMGW approximation, the screening is im- 
plicit in the self-consistency of charge densities in the 
HF product state. For weak and intermediate disorder, 
5min is quantitatively similar to the DFSGW case; how- 
ever, for strong disorder gmin < 1 for all nonzero U, in 
contrast to the DFSGW case. The quantitative similar- 
ity between the different (?min curves in Fig. [Ijc), unlike 
Fig. [Tfb), suggests that the self-consistent HF solutions 
screen the disorder potential more completely than the 
variational parameter Emm in the DFSGW approxima- 
tion. 

We now move to a discussion of the quality of the vari- 
ational solutions. In Fig. [21 we show the differences 
between the exact and variational energies. For weak 
and intermediate disorder (5-Egwf is smaller than (JE'uhf, 
suggesting that GWFs are better than UHF. Both the 
DFSGW and PMGW variational states have similar SE 
(recall that PMGW(g, e) is numerically equivalent to 
PMGW here), and detailed comparisons of the varia- 
tional states (not shown) find little difference in their 
predictions for the physical observables defined at the 
beginning of Sec. III. 

For strong disorder the situation is a little different: 
the value for (5£'gwf is different for the three variational 
states, with DFSGW having the smallest SE. However, 
for most values of U, the UHF approximation has a lower 
energy than any of the GWF approximations. 

In Figs [3] to [5] we compare charge densities and spin 
correlations for variational and exact calculations. For 
clarity, GWF results are shown only for the DFSGW 
approximation. Furthermore, we show the charge density 
as a function of U for only one representative site in the 
cluster, and we show the spin correlations for only a single 
representative pair of sites. 

Figure [31 shows the results for weak disorder, W/t = 
4/3. Here, W ^ D, and the bandwidth D is therefore the 
relevant energy scale for the crossover between the weakly 
and strongly-interacting limits as a function of U (similar 
to the disorder- free case''*). The figure shows that both 
the UHF and DFSGW do a good job of reproducing the 
local charge densities for all U, but that the DFSGW is 
significantly better at reproducing the spin-correlations. 
This is particularly true at large U where local moments 
have formed (i.e. all sites are singly occupied) : within the 
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FIG. 2; (Colour online) Difference between variational and ex- 
act energies for (a) weak, (b) intermediate, and (c) strong dis- 
order. Results are shown for DFSGW, PMGW(g, e), PMGW, 
and UHF. The UHF results are shown for step size AU = t. 
Note that for weak and intermediate disorder, the PMGW 
and PMGW((?, e) curves coincide. 
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FIG. 3: (Colour online) Comparison of (a) the charge density 
for the sixth site, ng, and (b) (Se ■ Si) for exact results (black 
dashed lines), DFSGW (solid red lines), and UHF(dotted blue 
lines). The UHF results are shown for step size Af/ = t. 
Resuhs are for W/t = 4/3. 



UHF approximation, local moment formation coincides 
with the onset of classical static magnetic moments, such 
that (Si ■ S,,+i) -1/4. In contrast, the DFSGW gives 
(Si • Si+i) —0.45 for large J7, which is the clean limit 
result for large C/— . 

Results for intermediate disorder are shown in Fig. 21 
and the conclusions are the same as for weak disorder: 
both the UHF and the DFSGW do a good job of re- 
producing the local charge density, but the UHF fails to 
reproduce the spin correlations. A closer examination re- 
veals that the UHF value for ng deviates from the exact 
fiQ near U = 1.3i, where local moments form at site 6. 

By contrast, results for large disorder in Fig. [5] show 
that while the DFSGW captures general trends, it does a 
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FIG. 4: (Colour online) As in Fig. H but for W/t = 4. 
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FIG. 5: (Colour online) As in Fig. El but for W/t = 12. 



poor job of reproducing the details of the charge densities 
and spin correlations as a function of U. Interestingly, 
the UHF does a remarkably good job of describing both 
the charge densities and spin correlations for U < 6t. We 
discuss the large-disorder limit in more detail below, but 
remark here that the success of UHF can be traced back 
to the fact that, for U <C W, most sites are empty or dou- 
bly occupied, with relatively few having moments. The 
UHF approximation correctly predicts that spin correla- 
tions between isolated moments and their near-neighbors 
vanish when the near-neighbor sites are empty or doubly 
occupied. 

Another measure of the quality of the various approxi- 
mations is the wavefunction overlap, shown in Figs. |S1 and 
III The results of these figures are essentially consistent 
with the results presented above: the DFSGW has a large 
overlap with the exact wave function (K^'exl^cwp)! > 
0.977) for weak and intermediate disorder, while for 
strong disorder the overlap is poor except in the small- 
and large-J7 limits. The UHF overlap is generally poor, 
except at small U. 

The reason for the difference in the quality of the GWF 
approximation at intermediate and large disorder is most 
easily understood by first considering the atomic limit 
(t = 0). When [/ = 0, the lowest-energy Ne/2 sites are 
doubly occupied while the remaining sites are empty. As 
U is increased, this arrangement persists until U is larger 
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FIG. 6: (Colour online) Magnitude of the overlap of the exact 
and DFSGW wave functions for the six-site disordered cluster. 



well for weak disorder {W ^ D) and fail for strong disor- 
der (VF D) where the physics approaches the atomic 
limit. The surprising result, evident in Fig. HI is that the 
moment formation is sufficiently homogeneous at inter- 
mediate disorder [W — D) to be well-represented by our 
simple variational wavefunction. 

In part, the success of the simple GWF approximation 
at intermediate disorder can be attributed to the screen- 
ing of the disorder potential by interactions (c.f. Fig. [T]). 
For small [/, correlation effects are minor and the wave- 
function is well represented by the original product state 
t/jps), while for [7 > D (where correlation effects are im- 
portant), Emin produccs a significant renormalization of 
the impurity potential. 
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FIG. 7: (Colour online) Magnitude of the overlap of the exact 
and UHF wave functions for the six-site disordered cluster. 
Data are shown for step size l\U = t 



than the energy difference between the highest-energy 
doubly-occupied site and the lowest-energy empty site, at 
which point one electron is transferred from the doubly- 
occupied site to the empty site. The process continues as 
U increases further, with electrons being promoted from 
doubly-occupied sites whenever the cost of double occu- 
pancy is greater than the cost of promoting an electron 
to the next available empty site. The process terminates 
when U > W, at which point all lattice sites are singly oc- 
cupied. The formation of local moments (singly-occupied 
sites) therefore happens inhomogeneously in the atomic 
limit. 

A nonzero t delocalizes electrons by an amount propor- 
tional to t/W, and thus makes the charge distribution 
more homogeneous. Given that the Gutzwiller projec- 
tion (which is responsible for generating local moments 
in I^gwf}) is chosen to be spatially uniform, it is un- 
surprising that the DFSGW approximation should work 



B. Ten-Site Cluster 



We now extend the work of the previous section to 
consider 10-site clusters with periodic boundary condi- 
tions. In this section, we consider results which are av- 
eraged over 10 randomly-generated complexions of the 
disorder potential. The site potentials are chosen to 
lie in the interval {—W/2,W/2) and we consider weak 
{W/t — 4/3), intermediate {W/t — 4), and strong dis- 
order {W/t = 12) cases as before. We show data for 
Ne ~ 10 electrons (half-filling) and — 6 electrons 
(near to quarter-filling). Given the large Hilbert space 
(63 500 states for Ng = 10), we use the Lanczos algo- 
rithm to find the exact ground states for comparison to 
the GWF results. 

Resuhs are shown for the DFSGW and PMGW ap- 
proximations, as well as for UHF. Results have not been 
shown for the PMGW(f/, e) approximation because it was 
found in Sec Hira to be identical to the PMGW approx- 
imation for weak and intermediate disorder, and because 
it is not significantly better than PMGW at large disor- 
der. 

Throughout this section, we show error bars for DF- 
SGW results. These error bars give the root mean 
square (rms) variation of Emm, ffmin, (SE), (Sn), and 
{dS ■ S') over 10 impurity configurations and (for (Sn) 
and {6S ■ S')) 10 sites. The bars are shown for every Ath 
data point and are of similar size in other GWF approx- 
imations. We emphasize that these error bars do not 
indicate the accuracy of the approximation (the curves 
themselves indicate this), but describe the site-to-site or 
sample-to-sample variation of the accuracy of the approx- 
imation. Thus, small error bars indicate that the approx- 
imation consistently over/underestimates a quantity by 
the same amount, while large error bars indicate that the 
quality of the approximation varies significantly from site 
to site. 
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FIG. 8: (Colour online) Variational parameters, (a) and (b), 
for the DFSGW and (c) PMGW approximations. Results 
are shown for ten-site disordered clusters averaged over ten 
complexions of disorder with A*'e — 10. Data are shown for 
step size AU = t. Error bars give the rms variation of Smin 
and Qmin over impurity configurations. 



FIG. 9: (Colour online) The disorder- averaged energy difi^er- 
ence per site between the variational and exact energies for 
(a) weak, (b) intermediate, and (c) strong disorder at half- 
filling {Ne = 10). The curves are DFSGW (solid black line), 
PMGW (red dashed line), and UHF (dotted blue line). Data 
are shown for step size AU = t. Error bars give the rms vari- 
ation of (SE) over impurity configurations, and are too small 
to see in (a) and (6). 



1. Ne = 10 electrons (half-filling) 



We begin by showing results for the disorder-averaged 
variational parameters as a function of U at half-filling 
(Fig. [8]) . The DFSGW curves are very similar to those 
shown in Fig. [1] for six sites. As before, differences be- 
tween the PMGW and DFSGW curves are most pro- 
nounced for large disorder and suggest that the PMHF 
approximation leads to a more complete screening of the 
disorder potential than docs DFSGW. 

The variational energies are shown in Fig. [5) For weak 
and intermediate disorder, the two GWF approximations 
give nearly identical values for ((5i?)GWF- The GWF en- 
ergies are lower than the UHF energies for small U, but 
larger than the UHF energies at large U. This is consis- 
tent with the results for larger clusters in the disorder-free 
case (c.f. Fig. 4 of Ref. 48), and is an indication that the 
simple GWF is unable to correctly reproduce the Mott 
transition at large U. For strong disorder, the Mott tran- 
sition occurs at U ^ W and this sets the energy scale at 
which the UHF becomes superior to the GWF. 

The error in the charge density (Sn), defined in 
Eq. ini), is shown in Fig. [TOl At half-filling, the re- 
sults are largely consistent with those found for six sites. 
Both the GWF and UHF approximations are good for 
weak and intermediate disorder, with the DFSGW ap- 
proximation producing an average error of less than 2.5%; 
however, all approximations work less well for strong dis- 
order. 

The spin correlations, defined in Eq. p^ . are shown in 
Fig. [TT] and are also similar to the six-site case. The UHF 
approximation does a poor job because it generates static 
local moments, while the GWF approximations work re- 
markably well for weak and intermediate disorder, hav- 




FIG. 10: (Colour online) A comparison of the mean abso- 
lute difi'erence between the exact and variational local charge 
densities for the ten-site disordered cluster averaged over ten 
complexions of disorder, with A^e = 10. Curves are for DF- 
SGW (solid black line), PMGW (red dashed line), and UHF 
(dotted blue line). Data are shown for step size AU = t. Er- 
ror bars give the rms variation of (Sn) over sites and impurity 
configurations. 



ing errors of order 1-2%. For strong disorder, the results 
are comparable for all approximations, when U < 7t. 
The GWF approximations are much better than UHF for 
U > 7. 
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FIG. 11: (Colour online) A comparison of the mean abso- 
lute difference between the exact and variational expectation 
values of the near- neighbour spin correlations for the ten-site 
disordered cluster averaged over ten complexions of disorder, 
with = 10. Curves are for DFSGW (solid black line), 
PMGW (red dashed line), and UHF (dotted blue line). Data 
are shown for step size AL'^ = t. Error bars give the rms 
variation of {5S ■ S') over sites and impurity configurations. 




FIG. 12: (Colour online) As in Fig. U but for A^^ = 6 elec- 
trons. 



2. Ne 



electrons 



Away from half-filling, strong correlations play a lesser 
role than at half-filling, and there is no Mott transition. 
The variational parameters for Ng = 6 are shown as a 
function of U in Fig. [121 The most noticeable difference 
with the half-filled case is that the variational parameters 
saturate at large U here. The energy scale at which satu- 
ration occurs appears to be D for weak and intermediate 
disorder, and W for strong disorder. A consequence of 
saturation is that the impurity potential is only partially 
screened at large U. Thus, while the half-filled ground 
state is the same for ordered and disordered models as 
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FIG. 13: (Colour online) As in Fig. U but for A^e = 6 elec- 
trons. 
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FIG. 14: (Colour online) As in Fig. [TOl but for Ae = 6 elec- 
trons. Error bars are larger here than for Ne — 10 because of 
the larger site-to-site variations in the accuracy of the GWF. 



U oo, the disorder potential remains relevant at all U 
for Ne = 6. This relevance is illustrated with a simple 
example: for W ^ t and U ^ oo, the ground state has 
the A'e sites with the lowest potentials Vi singly-occupied 
and the remaining N — empty. Since there is no Mott 
transition for A^g = 6, we expect the GWFs to be valid 
over a larger range of U than at half-filling. This is sup- 
ported by Fig. [131 where the GWF energies are generally 
lower than the UHF energies for weak and intermedi- 
ate disorder. As before, however, the GWF approxima- 
tions are worse at large disorder. It is interesting to note 
that ((5£')gwf increases roughly linearly with U at large 
U . Similar behavior was found previously in the ordered 
case,^* and is consistent with the fact that gmin is larger 
for Ne = Q than it is for — 10 (i.e. that the pro- 
jection of doubly-occupied states is less at A'g = 6 than 

Ne = 10). 

The charge densities for Ne = 6 (Fig. are similar to 
the half-filled case in the sense that the DFSGW charge 
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FIG. 15: (Colour online) As in Fig. [TT] but for iVe = 6 elec- 
trons. 



densities are comparable to UHF— at weak and interme- 
diate disorder, but are less accurate than UHF at strong 
disorder. Similarly to the Ne = 10 case, the average spin 
correlations (Fig. [15]) are reasonably well reproduced at 
weak and intermediate disorder. 



IV. DISCUSSION 




<5S.S'> 



FIG. 16: (Colour online) A comparison of the mean absolute 
difference between the exact and the variational energies (a), 
local charge densities (b), and near-neighbour spin correla- 
tions (c). The results are shown for ten-site clusters averaged 
over ten complexions of disorder (red solid lines for PMGW 
and blue dot-dashed lines for UHF), and six-site clusters av- 
eraged over fifteen configurations of disorder (diamonds for 
PMGW and blue dotted lines for UHF). Data are shown for 
step size Ai7 = t for ten-site clusters and AI7 = Q.lt for 
six-site clusters. 



The GWFs used in this work make two key simplify- 
ing assumptions: first, that the correlation physics is ap- 
proximately local and can be represented by a Gutzwiller 
partially-projected wavefunction and, second, that the 
spatial inhomogeneity associated with the disorder po- 
tential can be incorporated entirely within the uncorre- 
lated product state. In this section, we discuss these two 
assumptions, keeping in mind that, while it is possible 
to treat a large number of variational parameters using 
statistical sampling methods,-^ many of the interesting 
questions relating to disordered systems require the abil- 
ity to study large system sizes, and that it is therefore 
desirable to restrict the number of variational parameters 
as much as possible. The results shown in this paper give 
some clues as to how this can be done. 

There are two distinct regimes in our calculations, dis- 
tinguished by the strength of the screened disorder po- 
tential, which in the DFSGW is W = Wjs^,^. For 
< D, W is sufficiently small that the GWFs used in 
this work give similar results for {SE) as existin g ca lcu- 
lations for the ordered case (i.e. Fig. 4 of Ref. 'Is). In 
Ref.[48|, it was shown that the simplest GWF could be im- 
proved by inclusion of a long-range Jastrow factor, and 
the similarity of our results to theirs has a similar im- 
plication. Interestingly, there is nothing in our results 
to suggest that adopting spatially-inhomogeneous varia- 
tional parameters would make a significant improvement. 
(Recall, Figs. 13] and m the accuracy with which the DF- 
SGW approximation reproduces the local charge density 
and spin correlations.) It appears as if, for W 



spatial inhomogeneity can be adequately incorporated in 
the product state provided that screening is included, ei- 
ther through a Hartree-Fock determination of the prod- 
uct state or through a variational parameter. This rep- 
resents a huge savings of computational effort since, if 
one makes a physically reasonable ansatz for the form 
of the Jastrow factor, the number of variational param- 
eters can be reduced to four: a screening parameter e, a 
projection parameter and a pair of parameters for the 
Jastrow factor.— 

For strong disorder, the situation is different since spa- 
tial fluctuations of the charge density are significant. It 
appears that, in this case, spatial inhomogeneity cannot 
be included solely in l^'ps), but must be included in the 
variational parameters as well. This was the approach 
taken in Ref. [s^ where the local variational parameters 
g and e were allowed to vary from site to site (but the 
nonlocal Jastrow factors were treated as spatially homo- 
geneous in order to reduce the computational workload) 
on lattices of TV ~ 100 sites. It is unclear that any simpli- 
fication of this approach is possible in the large-disorder 
limit. 

Finally, we comment on finite size effects. We compare 
(Jn), and ((5S • S'), for six and ten site clusters at 
half filling and intermediate disorder in Fig. [1^] The 
GWF energy differences increase with system size and, 
in the weak disorder case, appear to scale towards the 
clean limit results shown in Fig. 4 of Ref. |4^ for large 
systems. We note that, while there are quantitative dif- 
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ferences , (SE) for the 10-site clusters has the same qual- 
itative U-dependence as for the larger clean limit sys- 
tems. More importantly for this work, (Sn) and (SS ■ S') 
are nearly independent of system size. This is consistent 
with earlier scaling results in the disorder-free case^^ and 
suggests that the high accuracy of the results in Fig. [11] 
(b) and (c) will extend to larger systems. Thus, ability to 
benchmark GWF trial functions is not inhibited by finite 
size effects. 



V. CONCLUSIONS 

We have studied the Anderson-Hubbard model, one 
of the simplest model Hamiltonians that can describe 
correlated electrons moving on a disordered lattice. We 
have compared the exact and variational ground states 
for disordered one dimensional chains up to a length of 
10 sites. The main focus of our calculations was on find- 
ing the simplest variational states that accurately include 
the effects of disorder. The quality of the approxima- 
tions was determined by comparing a number of physical 
quantities — the ground state energy, local charge densi- 
ties, and near-neighbour spin correlations — to the exact 
results. Because it has been, until recently, the standard 
tool for studying disordered systems, we also compared 
our variational states to the unrestricted Hartree-Fock 
approximat ion . 

Our main conclusion is that for weak and intermediate 
disorder, an accurate description of the ground state can 
be obtained by taking a spatially inhomogeneous product 
state, and spatially homogeneous variational parameters. 
This represents a significant reduction of computational 
workload over the more general case where both the prod- 
uct state and variational parameters must be inhomoge- 



neous. 



APPENDIX A: TWO-SITE CLUSTER 

We compare the different product states for a two-site 
system, where an analytical solution can be found. For 
this cluster we choose the on-site energy of sites one and 
two to be V and respectively, and take — 2 (corre- 
sponding to half- filling) . Then 



'ps/ 



1 



(Al) 



y^l + +a2 + (a|aJ2 
x(| n)|0) - atlDI T) + aj T>U> + atailO)! Ti)), 

where N„~^^^ ^ are normalized eigenvectors of the 
matrix form of Hbi for spin cr. 



with V = V/e, 



-t Un2-a 



(A2) 



Oct = -^{V + U{ni^a -n2- 



+ - n2-a) + V'Y + 4^2}, (A3) 

and the many-particle states are defined as | ti)|0) = 
cItCu|0)> I t)l i) =cj^4il0), etc. 

Using Eqs. ^ and \k\\ . \^gwf) is written as 



\'^GWf) = . . — (5l ti)|0) - a^l i)| T) +aj T)l i) ^ ga^am \\)) , (A4) 

y'l + + (a^aiY 



and 

.g2[(l + a2a2)[/ + 2V\ - 2tg{a^ + ax)(aTai + 1) + ^("t + "i) 
5'(l + (aTai)') + a^+aJ 

r 



^GWF = '■ 1277^:77—^2^—^—2 • (A5) 



Numerical results for this cluster show that the DF- 
SGW, PMGW(5,e), and PMGW energies are the same, 
and we therefore only show the results for PMGW in 
Fig. [T71 In this figure, we compare the magnitude of the 
PMGW energies with the UHF and exact energies. We 
also compare the overlaps of the PMGW and UHF wave- 
functions with the exact wavefunctions. The overlaps are 
defined as |(^ox|*var)|, where (^'ox| is the complex con- 
jugate of the exact ground state eigenvector of H. The 



reason to study both ground state energies and overlaps 
is that a lower variational energy doesn't necessarily im- 
ply a better variational wave function. 

In general, UHF energies are lower than PMHF ener- 
gies (not shown), suggesting that the UHF approxima- 
tion is better than the PMHF one. However, after the 
projection the situation is reversed. In Fig. [T7| one sees 
that PMGW energies are exact, and that the UHF ener- 
gies are higher than the PMGW energies for U/t > 2, 
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where magnetic moments develop in the UHF states. 
Similarly, \{'i'ex\'^ pmgw)\ = 1 while |(«'ea.|4'uHF)| de- 
creases rapidly as moments develop. This is consistent 
with the ordered case,^® where the GWF is exact for the 
two-site problem, but is approximate for larger system 
sizes. 
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FIG. 17: (Colour online) Variational solutions for the two-site 
Anderson-Hubbard model with V = 2.0t. (a) a comparison 
of the PMGW variational energy with the exact and UHF 
ground-state energies. The exact energies coincide with the 
PMGW energies, (b) the wavefunction overlap |{\E'cx|*I'var)|. 
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